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Abstract 

Gravitational waves (GWs) are one of the key signatures of cosmic strings. If GWs from 
cosmic strings are detected in future experiments, not only their existence can be conhrmed but 
also their properties might be probed. In this paper, we study the determination of cosmic string 
parameters through direct detection of GW signatures in future ground-based GW experiments. 
We consider two types of GWs, bursts and the stochastic GW background, which provide us with 
different information about cosmic string properties. Performing the Fisher matrix calculation 
on the cosmic string parameters, such as parameters governing the string tension G^i and initial 
loop size a and the reconnection probability p, we find that the two different types of GW can 
break degeneracies in some of these parameters and provide better constraints than those from 
each measurement. 
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1 Introduction 



Cosmic strings are one dimensional topological defects which can be formed at a phase transition in 
the early Universe [1] (as a review, see [2]). Strings in superstring theory which have cosmological 
length, so-called cosmic super strings, can also appear after the stringy model of inflation and behave 
like cosmic strings [31 IH |5]. They form the highly complicated string network, including infinite 
strings, which stretch across Hubble horizon, and closed loops, whose size is much smaller than the 
Hubble scale, and leave cosmological effects in many ways through their nonlinear evolution. They 
therefore have attracted strong attention since the possibility of their existence was pointed out. 
Considerable research concerning their observational signals and their detectability have been done 
so far. If their signals can be observed precisely enough, not only the existence of cosmic strings 
will be confirmed but also their properties might be studied. Since we can probe physics beyond the 
Standard Model of particle physics such as grand unified theory (GUT) or superstring theory through 
the study of cosmic strings, it is very interesting and meaningful to examine how the properties of 
cosmic strings can be constrained through future experiments and observations. 

Among important signals of cosmic strings are gravitational waves (GWs) EllZl El El [101 IHl [121 
[131 [HI [151 [IS] • The main source of GWs in the string network is cusps on loops, [j A cusp is a highly 
Lorentz boosted region on a loop which appears 0{1) times in an oscillation period of the loop. 
Beamed GW bursts are emitted from cusps and can be detected directly as strong but infrequent 
bursts as well as in the form of a stochastic GW background, which consists of many small bursts 
overlapping each other [101 [H] . The rate of GW bursts and the spectrum of the GW background 
depend on the parameters which characterize the string network. Conversely, we can constrain the 
parameters from the fact that GWs from cosmic strings have never been detected, or even we can 
determine the values of them through future observations if GWs from cosmic strings are detected. 

So far, constraints on cosmic string parameters have been imposed by LIGO from non-detection 
of either bursts or GW background [101 [20]. Currently LIGO is undergoing a major upgrade (Ad- 
vanced LIGO) [21] and will increase the detector sensitivity by more than an order of magnitude. 
Furthermore, several new ground-based detectors would be built along the same time line as LIGO, 
such as LCGT in Japan [22], Advanced Virgo in Italy [23], and AIGO in Australia [21]. These 
additional detectors form a world-wide network of gravitational- wave observatories, which can be a 
more powerful tool to search for signatures of cosmic strings. 

Cosmic strings can be characterized by the following parameters. The most important one is 
the tension /i, the energy stored per unit length in a cosmic string, which is often written in the 
form of its product and Newton constant, Gfi. For field theoretic cosmic strings, fi is roughly the 
square of the energy scale at the phase transition which leads to the appearance of cosmic strings. 
If they appear at the spontaneous breaking of the symmetry of GUT, expected value of G/i is of 
O{10~^) or (9(10^^). For cosmic superstrings, /i is proportional to the square of the string scale, but 
it also depends on other parameters such as the warp factor of the extra dimension where the cosmic 
superstring is located. Therefore, the tension of cosmic superstring can take a broad range of values. 
Since Gfi determines not only the typical amplitude of GWs emitted from loops but also the lifetime 

^ GWs from kinks may dominate that from cusps in the case of cosmic superstring |17[ 118) . However their con- 
tribution strongly depends on the fraction of loops with junctions. So we do not consider their contribution in this 
paper. 
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of loops, it affects the amplitude of GW bursts and the GW background as well as the rate of bursts 
and the spectral shape of the GW background. 

The second one is the loop size a. It is well known that infinite strings reach the scaling regime, 
where the curvature radius and the interval of strings are comparable to the Hubble radius. They 
have to lose their length by releasing loops continuously in order to maintain the scaling. It is 
considered that the typical size of initial loops also obey the scaling law and it is often written in 
the form of at. Although many analytic or numerical studies have been conducted to determine the 
value ofa|25l[26l|271[28l|29lliol[3ll|33[33l^ inferred values of a range from 0.1 to some 

powers of GyU, so its magnitude remains controversial. We treat it as a free parameter in this paper. 
This parameter a also affects the rate of GW bursts and the spectrum of the GW background, since 
lifetime of loops and frequency of GWs emitted from loops depend on the loop size. 

The third one is the reconnection probability p. It is well known that, in almost all models, 
field theoretic cosmic strings necessarily reconnect when they collide. On the other hand, cosmic 
superstrings can have the reconnection probability much smaller than 1, because a collision between 
them is a quantum process and they can miss each other in the extra dimension [36| W(\ [38] . The 
reduced reconnection probability leads to an inefficient loss of length of infinite strings and eventually 
causes an enhancement of their density [SHI SO]- Since the number of loops accordingly increases, the 
amplitude of GW background and the burst rate are enhanced as p decreases. 

In this paper, we calculate the spectrum of the GW background and the rate of "rare bursts", 
which are isolated and whose amplitudes exceed that of the background, and study their dependen- 
cies on above three parameters. Then we find the parameter region which is excluded by current 
experiments or can be searched by future experiments. Finally, assuming that GWs will be detected 
by ground-based GW detectors, we use the Fisher matrix formalism in order to investigate the degree 
to which the cosmic string parameters are constrained by future experiments. It is notable that the 
rare bursts and the stochastic GW background have different information and constraints from them 
break the parameter degeneracies each other. 

This paper is constructed as follows. In the next section, we show the model of the string network 
assumed in this paper. In section 3, we describe the formalism for calculation of the burst rate and the 
GW background spectrum and show some examples assuming specific parameter values. In section 4, 
we find the parameter region which can be probed by upcoming GW experiments such as Advanced 
LIGO and the world wide GW network. Furthermore, we calculate the Fisher information matrix and 
predict the combined constraints from burst detection and GW background measurements, choosing 
parameter values where both are accessible by future experiments. The last section is devoted to the 
summary. 

2 Analytic model of cosmic string network 
2.1 Infinite strings 

We adopt the model in [IHl E], which is based on the velocity-dependent one-scale model [12]. The 
network of infinite strings can be considered as a random walk, so their total length L in volume V 
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Figure 1: The values of 7 plotted as a function of p for the radiation and matter-dominated era. 



can be written as 



(1) 



where ^ is the correlation length, which corresponds to the typical curvature radius and interval of 
infinite strings. The equations for 7 = ^/t and the root mean velocity of infinite strings v are given 
by 

i^_l+.+|^; + .„^ (2) 

7 at 27 



dv 
It 



^ = (1 



.k{v) 
z/7 



2v), 



(3) 



where k{v) = i+sS ' ^ Hubble parameter, and the scale factor a is parameterized as 

a{t) oc f^. The constant parameter c represents the efficiency of loop formation and we set it to 
be c = 0.23 according to Ref. [H]. In the scaling regime, 7 and v become constant. We can get 
their asymptotic values by setting d'y/dt and dv/dt to be 0. Figure [1] shows the value of 7 in the 
radiation-dominated era {u = 1/2) and the matter-dominated era {u = 2/3) as a function of the 
reconnection probability p. This figure shows that 7 is proportional to p and p^^^ in the radiation 
and matter-dominated era, respectively. Hereafter, we denote the values of 7 in the radiation and 
matter-dominated era as 7,. and 7^- Since the effect of the time dependence of 7 around the matter- 
radiation equality is small for the parameters we set below, we approximate 7 as a step function. 



7(2;) 



7r 



:z > z, 



eq 



7™ 1 ^ ^ ^eq 



(4) 



where Zgq is the redshift at the matter-radiation equality. 
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2.2 Loops 

In order to maintain the scaling, infinite strings have to continuously release their length in the form 
of loops. The released length in a Hubble volume per Hubble time is comparable to the length of 
infinite strings in a Hubble volume. Since the length of a loop formed at time t is given by at, the 
number density of loops produced between time t and t + dt is 

n{t)dt = (5) 

The number density of loops is diluted proportional to by cosmic expansion. Therefore, the 
number density of loops formed between ti and ti + dti at time t is 

nit,U)dU = ^J^y . (6) 



a'-fHf \ a{t) 

After loop formation, it continues to shrink by emitting its energy as GWs. The energy spectrum of 
GWs from a loop of circumference / per unit time is given by 

^ ~ G^^H-'/'f-'/' (7) 

with a low frequency cutoff at / ~ l^^. The total energy emission rate is 

E = TGfi^ (8) 

where F is set to be 50 in this paper. Then, the length of a loop formed at ti is 

lit,ti) = ati-TGfi{t-ti), (9) 

at time t. Lifetime of a loop formed at ti is given by f^ti. If a < TGfi they are short-lived, that 
is, they decay within a Hubble time. If a > TGfi, loops are long-lived, that is, they live longer than 
a Hubble time. In this case, loops have a wide range of length from at to TGfit, which corresponds 
to loops just formed and expiring at some time. The most numerous loops are those of length 
comparable to ~ TGfit. 



3 Gravitational waves from cosmic string loops 

In this section, we describe the calculation procedure for the burst rate and the spectrum of the GW 
background, following the formalism in Ref. [TT]. Using the result, we show parameter dependence 
of the burst rate and the background spectrum for several parameter sets. 

3.1 Formalism 

The linearly polarized waveform of a GW burst emitted in a direction n by a loop with circumference 
/ at redshift z is expressed as 

V(i,n)= [ dfhif,z,l)e-^^'f'eM X Qi^-^c - cos9Uf,z,l)) X Q{1 - 9Uf,z,l)), (10) 
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where ric is the direction of the center of the burst, which coincides with that of the derivative of the 
right or left moving mode of the loop, and 6^ is the beaming angle of the GW burst which is given 
by 

9Uf,z,l) = {{l + z)fl)-'/'. (11) 

The first Heviside step function 6 in Eq. (11 01) reflects the fact that the burst is emitted into a 
limited angle, and the second one means that it has a low frequency cutoff at / < /"^ when it is 
emitted. The polarization tensor (for plus polarization) is expressed as e^^ = l^m^ — luin^, where 
Ifj, = (0, 1), = (0, m) and 1 and m are unit vectors orthogonal to n and each other. 

A GW burst from a loop is most efficiently generated at frequency comparable to and the 
amplitude of higher frequency modes decrease in proportion to /~^/^[^. The Fourier transform of the 
GW amplitude is given byH 

47r(12)^/3 Q^i 
' ^ {2ny/^{3T{l/3))mi + z)fiy/^r{z)f ^ ' {{1 + z)fiy/^r{z)f' ^ ^ 

where 



Note that one can confirm that this frequency leads to the energy spectrum of GWs from a loop 
given in Eq. ([7]) (see, for example, Ref. [13]). 

The number of GWs coming to the Earth per unit time, emitted at redshift between z and z + dz 
by loops formed between ti and ti + dti i^ 

-dzdti = -^OrrXf, Z, lf——^^—-—^{t{z),U)dU^dz X 9(1 - e„(/, Z, I)). (14) 



dzdti 4:71 (1 + z)l{t{z),ti) dti dz 

Here, -^Omifyzjy is the factor which represents the fraction of GW bursts beamed towards the 
Earth. Cusp formation is expected to occur 0{1) times in an oscillation period, which is characterized 
by parameter c. We set it to be 1 in this paper. The factor ^dz is the volume between z and z + dz 
at time t{z) and given by 

dV Ana'izyjz) 
dz^^'~ H{z){l + z)- 

Since the observables are the amplitude and the frequency of GW bursts, we need a prediction 
of the burst rate expressed in terms of the given amplitude and frequency. Using Eqs. iQ and f|T2|) . 
we can rewrite Eq. f|T^ to express the number of GWs coming per unit time which were emitted at 
redshift z and which have frequency / and amplitude h at the present time 

'^^if,h,z) = ^el{f,z,i)-^^,^-— ( ^X^Q{i-euf,zJ)), (16) 



dzdh ' ' 471"^'' {1 + z)h'-f'^atj a + TGfi \a{t) J dz 



^ Strictly speaking, Eq. (jl2p is valid only for modes whose frequency is much higher than when the GW 

is emitted and the amplitude of low frequency modes, / ~ depends on the detail of the oscillation of the loop. 
However, it is known from numerical studies that it is a good approximation to apply the power law as in Eq. (|12|) to 
low frequency modes for small loops [48]. 

^ Note that this definition h is the same as Ref. [T3] and different from Ref. [11] by a factor of /. 
We include another factor i , which is omitted in Ref. [TT] . This leads to more accurate estimation of burst rate. 
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where I and ti can be given as functions of /i, z, and /, 



Kf,h,z)=i,!^{l + z)'/'fA'^\ (17) 



U{f, h, z) 



2.68G/i' ' ' ) 
l{f, h,z) + TGfit{z) 



a + TGfx 

By integrating Eq. (fT6|) in terms of z, we get the rate of GWs for the given frequency and amphtude, 

^= r dz^ (19) 

dh Jq dhdz 

Bursts overlapping each other form a stochastic background. We adopt the criterion in Ref. [T3] . 
which counts bursts coming to the Earth with a time interval shorter than the oscillation period of 
themselves as a component of the GW background. Such bursts have amplitude smaller than /;,*, 
which is determined for a given frequency as 

K dh Jq dzdh 

The amplitude of a stochastic GW background is commonly expressed by flcwif) = {dpQ-^/dln f)/ 
where pcw is the energy density of the GWs and pcr is the critical density of the Universe. Then 
the spectral amplitude of the GW background is given by summing up all contributions from small 
bursts, 

27r2 r''* ,_,dR 



ficw(/) = ^/»| dhh^^. (21) 

3.2 Parameter dependence of the burst rate and the GW background 
spectrum 

Here, we show how the burst rate dR/dh as a function of amplitude h depends on parameters Gfi, 
a, and p. In Figure [21 we plot dR/dlnh versus h for Gfi = 10""^, a = 10~^^, p = 1, which is evaluated 
at / = 220Hz, the best-sensitivity frequency of Advanced LIGO. 

We set cosmological parameters, to be the WMAP 7-year mean values |19]: the dark energy 
density divided by the critical density Q\ = 0.728, the dark matter density divided by the critical 
density = 0.272, the current CMB temperature Tq = 2.725K, and the Hubble constant Hq = 
70.4km/s/Mpc. We see the natural tendency that stronger bursts have a lower rate. In this figure, we 
also show the ranges of h which correspond to bursts observed as rare bursts or the GW background 
by Advanced LIGO. GW bursts whose amplitude are larger than the detector sensitivity, which 
corresponds to fh = 3.4 x 10"^'^ at / = 220Hz for Advanced LIGO, are observed as rare bursts. In 
contrast, GW bursts whose rate is larger than their frequency, / ~ 220Hz at Advanced LIGO, are 
measured as the GW background, [f] 



^Note that bursts in the middle amplitude range between the rare burst and GW background regions may be 
detected as unresolved sources in the GW background [H]. Non-gaussian measurements of the GW background may 
be useful to characterize their contributions [151 [13 [17] . 
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Figure 2: The rate of GW burst dR/dlnh as a function of fh, for Gfi = 10"'', a = 10"^^, p = 1, 
and / = 220Hz, represented by the black thick hne. Bursts in the left region (orange) form the GW 
background and those in the right region (blue) are observed as an isolated burst. The red, green 
and purple dotted line represent the contribution from 1 < 2; < 10, 0.1 < 2; < 1 and 10"^ < < 0.1, 
respectively. 



GW bursts in each amplitude bin consist of bursts from different redshifts. In order to illustrate 
which redshift mainly contributes to bursts at a given amplitude, we also plot the contributions to 
dR/dln h from different redshift ranges in Figure 121 The red, green, and purple dotted line correspond 
to bursts from 1 < 2; < 10, 10~^ < 2 < 1, and 10"^ < < 10"^, respectively. This indicates that, in 
this parameter set, bursts detectable by Advanced LIGO come from redshift lower than z ~ 10~^, 
and the GW background consists of bursts emitted at redshifts higher than 2 ~ 1. 

show dependencies of the burst rate on parameters. 



Figures 3(a), 3(b), and 3(c 



Gjj, a, and 

p, respectively. The parametric dependence of the burst rate is the key in studying constraints 
on cosmic string parameters by GW direct detection experiments, since it determines the direction 
of the parameter degeneracy. The details are investigated in Appendix |X1 Here we only give some 
short explanations. The dependence on p is simplest because a small value of p simply enhances 
the rate through the factor 7"^ in Eq. ( fT6|) . A large value of G/i basically leads to a larger burst 
rate because of the enhancement of the amplitude of each burst. However the actual dependence 
is more nontrivial since the variation of Gfi also changes the typical lifetime of loops, which leads 
to difference in their density. Variation of a also affects the lifetime of loops, as well as the initial 
number of loops. Therefore the dependence on a is also not easy to simplify. 

We also show the spectrum of the GW background, flew, for different parameter sets in Fig. 



4(a) and Fig. 4(b) We also explain details of the parameter dependence of I^gw in Appendix [Bl 
Physically, ^Iqw at each frequency reflects the energy ratio of expiring loops to the total energy of the 
Universe at the redshift corresponding to min{a, FG/i} xt ~ f~^{l+z)~^ for the radiation-dominated 
Universe, or, for low frequency modes, it reflects the energy of GWs emitted recently. 
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fh 

(c) For Gfi = 10-^ a = IQ-^^, f = 220Hz, and various p. 
Figure 3: The burst rate dR/dlnh in terms of fh for various parameter sets. 
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(a) For G/i = 10 ^, p = 1, and various a. (b) For a = 10^^^, P = 1, and various G^i. 



Figure 4: The background spectrum ^qw versus present frequency / for various parameter sets. 

For these reasons, rare bursts and the background spectrum carry information of cosmic strings 
at different redshifts, and have different dependence on cosmic string parameters. This is one of the 
reason they provide different directions of parameter degeneracy in parameter constraints by burst 
detection and the GW background measurements, as we shall present in the next section. 



4 Constraints on cosmic string parameters from future GW 
experiments 

In this section, we investigate how accurately the cosmic string parameters can be determined if 
GWs from cosmic strings are detected by future experiments. We consider the case where both the 
GW bursts and the GW background are detected and show their constraints complement each other. 
Before that, we compute the sensitivities of GW detectors and show the accessible parameter space 
by current and future experiments. 

4.1 Sensitivity for GW detection 

In this paper, we assume the world wide GW network, consists of Advanced LIGO pairs, LCGT, 
VIRGO, and AIGO. Additional detectors improve angular resolution and enable precise determi- 
nation of intrinsic amplitude and polarization of the incoming GWs. Also, since these world wide 
GW detectors are directed to different regions of the sky, they provide almost all sky coverage and 
increase the number of detection events, while each GW detector has limited sky coverage. The 
detector network also improves sensitivity of the GW background search with the cross-correlation 
analysis. 
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4.1.1 Burst 



As shown in Eq. (fT2|) . the burst signal from a cosmic string cusp would be linearly polarized and 
have the frequency dependence of f~^^^ with cutoffs at low and high frequencies. We assume that 
the spectrum has the form of 

= Ar'/'e{u - /)e(/ - fi), (22) 

where the amplitude A can be read from Eq. f lT^ . The low frequency cutoff corresponds to the size 
of the loop, whose scale is typically cosmological. So, usually, the cutoff frequency is much lower 
than the lower frequency limit of direct detection experiments. The high frequency cutoff depends 
on the viewing angle as fh ~ 2/{9l^^L), where 6'obs is defined as the angular separation between the 
observer's line of sight and the beam direction of the GW emission. 

The total output of the detector is written as a combination of GW signal h(t) = F^h^{t) and 
detector noise n(t), s{t) = h(t) + n{t), where F"*" describes the detector response to plus polarized 
GWs. This function depends on the sky location of the GW source, the polarization angle and the 
configuration of the detectors, and can be interpreted as the sky area covered by the experiments. For 
a single detector, one may use the all sky-averaged value for orthogonal arm detectors, F+ ~ 1/ y/E. 
For next generation detector network (Advanced LIGO, LCGT, VIRGO, and AIGO), we assume 
that the GW detector network has 100% visibility over the whole sky so that -F+ ~ 1. 

The search for GW bursts signals from cosmic strings is usually performed via matched filtering 
[50| [T9] . The template for cosmic string bursts is usually taken as 



T 



{f) = r'''Q{h-f){f-fi). (23) 



The template is normalized by dividing by a = \/Jt\t), where the inner product is defined as 

[x\y)^A^j^df-^^y^. (24) 

Thus the normalized template is f = r/o", which satisfies {f\f) = 1. The noise spectral density Sn{f) 
is defined by (n(/)*n(/)) = 5'„(/)5(/ — f')/2, where (■ ■ ■ ) denotes the ensemble average. The noise 
spectral density for Advanced LIGO is given by [5T] 



Snif) = 10- 



-49 



Hz-\ (25) 



where x = f/ (215Hz), and the noise for current LIGO is given by |50j . 



28 / T nriTT.. \ 4 

S^{f) = 1.09 X 10-"^^ ( ) + 1.44 X 10-^5 



/100Hz\' 

1,28 X 10-« I 1 + ('^5^5^') 'jHz-'. (26) 



For other GW network observatories, we assume all detectors have the same sensitivity as Advanced 
LIGO. 
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The signal to noise ratio p is given by the product of the signal and the normalized template as 
p = (s|f). This is equivalent to calculate 




The low frequency cutoff is determined by the detector's limitation, which is taken to be // = lOHz 
for Advanced LIGO and fi = 40Hz for current LIGO. The high frequency cutoff is different for 
each burst, depending on the viewing angle. However, in most cases of detection, it is distributed 
around the most sensitive frequency of the detector [SI]. Here, we take fh = 220Hz for Advanced 
LIGO and fh = 150Hz for current LIGO. Taking the detection threshold as p > 4 [TU], we find 
the future GW detector network (Advanced LIGO sensitivity with full sky coverage -F+ ~ 1) can 
detect GW bursts whose amplitudes are larger than A ~ 2.1 x 10~^^ s~^/'^, which corresponds to 
fh ~ 3.4 X 10-23 at / = 220Hz. For current LIGO (LIGO sensitivity with F+ ~ 1/^5), the detection 
limit is A ~ 9.1 x 10"^^ s"^/^ which corresponds to fh ~ 1.7 x 10"^^ at / = 150Hz. 



4.1.2 stochastic GW background 

The GW background is searched by correlating output signals of two or multiple detectors. One may 
define the cross correlation signal between two detectors labeled by / and J as [52] 

rT/2 j-T/2 

S= dt dt'si{t)sj{t')Q{t,t'), (28) 

where T is the observation time and Q{t, t') is a filter function. Using the fact that noises of different 
detectors have no correlation each other, {si{t)sj{t')) ^ {hj{t)hj{t')), and transforming to Fourier 
space, the mean value of the signal can be expressed as 

/•oo roo 

p^{S)= df df'Srif - f'm{f)h.j{f'))Q{f'), (29) 



where the tilde denotes Fourier-transformed quantities and Sxif—f) = '^'^^^^ ~ ^^^{'^ fT) / n f . 

The response of the detector is given using as 

hiif) = E / ^^^^(/' ^)e-'^'^^-'''F^U\ (30) 

A ^ 

where A runs for both plus (+) and cross (x) polarization and x/ denotes the position of the detector. 
Using the relation between the Fourier amplitudes hx{f, f2) and flaw, 

(/^l(/,^2)/^v(/^^^')) = ^,5'{n,n')l5xx'S{f-f')\f\-'ncwm, (si) 

the cross correlation signal is given by 

P = / df\friij{f)^Gw{f)Q{f), (32) 
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where we define the overlap reduction function as 



7/j(/) = ^ 
8n 



(33) 



We calculate the overlap reduction function following the procedure given in Ref. whose Table 
2 or 3 provides the relative positions of future ground-based GW detectors. In the weak-signal 
assumption, the variance of the correlation signal is 



a 



is')- 

T/2 



T/2 



T/2 

dt I dt'{sjit)sjit)sjit')sjit'))Qit)Qit') 

T/2 



(34) 
(35) 



Using {sj{t)sj{t)si(t')sj(t')) ^ {nj(t)nj(t')) {nj(t)nj{t')) and transforming to Fourier space, this can 
be expressed in terms of the noise spectral density Sn{f) as 



0" ~ — 



#^n,/(|/|)^n,j(|/|)|Q(/)P. 



(36) 



The signal to noise ratio is defined hy p = n/a. Choosing the optimal function to maximize the 



signal to noise ratio, which is Q{f) cc \fll3s^^'l{\^\)s'^\^{\f\) ' obtain 



PiJ 



3Hl 
IOtt^ 



df 



uo 



1/2 



(37) 



For a network of detectors, we can make N[N — l)/2 independent correlation signals. One can 
calculate signal to noise ratio as (see Sec. V-C of Ref. [52] ) 



' N N 
.1=1 J<l 



1/2 



(3J 



Advanced LIGO detector pair would be able to detect the GW background with p > 4 if flaw > 
5.1 X 10~^ (for a flat spectrum) with 3- year observation. The multiple detector network would 
reach ^gw ~ 3.6 x 10~^, where we take the range of integration from lOHz to 3000Hz. The cross- 
correlation analysis with two-year run of the current LIGO detector pair has placed an upper limit 
Oqw < 7.2 X 10"^ which was performed for the frequency band 41.5Hz-169.25Hz. 



4.2 Accessible parameter space of cosmic string search 

Here, we find the parameter space excluded by current GW experiments and cosmological constraints. 



and that accessible by future GW experiments. In Fig 5(a), 5(b), and 5(c), we show the parameter 
regions which are excluded or can be probed by GW experiments in the a — Gfi plane for p = 1, 
10~^, and 10~^, respectively. Current constraints are provided by pulsar timing and current LIGO 
experiments. We also show the cosmological constraints from CMB and BBN, which comes from 
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current pulsar 

,-11^ LIG0S5GWB 



LIG0S5 burst 
AdvLIGO GWB 
AdvLIGO burst 




(c) p = 10-2 



Figure 5: The parameter regions excluded by current experiments and cosmological constraints, as 
well as regions which can be probed by Advanced LIGO, for different values of p. The colored regions 
are excluded by current experiments. Advanced LIGO can probe the region above the yellow solid 
line with GW background search and that above the yellow dotted line with burst search. 
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the fact that the energy density of the GW background, which has the same effect as extra neutrino 
species on CMB and BBN, must be small at the last scattering epoch and BEN, so as not to distort 
the fluctuation of CMB and change abundance of various nuclei. For interferometers, we show both 
parameter regions accessible by burst and GW background search with GW detector network consists 
of Advanced LIGO and other next generation detectors. For burst detection, we define the detection 
criterion to be whether the bursts whose amplitude is equal to the detector's best sensitivity come 
with a rate higher than lyr~^. Following the discussion in the previous section, we take the best 
sensitivity of current LIGO as //i = 1.7 x 10^^^ at / = 150Hz, and the best sensitivity of Advanced 
LIGO with GW detector network as fh = 3.4 x lO'^s at / = 220Hz. For the GW background 
search, we assume it is detectable if the amplitude is higher than Qqw = 7.2 x 10~^ for LIGO, and 
Qgw = 3.6 X 10-9 for Advanced LIGO. 

For the upper limit of ^Iqw from current pulsar timing experiments, we take 1.9 x 10"® at 
/ = 3.2 X IQ-^Hz [55]. The CMB provides the constraint that / r2Gw(/)'^(ln /) must be less than 
1.4 X IQ-^ at the last scattering [56]. The BBN constraint is / flGw{f)d{^^ f) < 1.6 x 10"^ at the epoch 
of BBN [131 EZ] • The lower limit of the integral is determined by the lowest frequency of the GWs 
emitted by largest and youngest loops at the time of CMB and BBN. The upper limit is the frequency 
of GWs emitted by the earliest loops when they expire, which is given by / ~ j^CrGfit(zi))~^, where 
z is the redshift at the last scattering or BBN and Zi is the redshift when such loops expire. We 
consider that such loops are formed at the end of the friction domination, when the temperature of 
the Universe is ~ y/Gfi. 

The reason why the current pulsar timing experiments constrain only the large a region is that, 
for a small value of a, loops cannot generate GWs at such low frequencies accessible by pulsar 
timing experiments, which is comparable to lyr"^. We see that although current constraints on 
the parameter region are rather severe especially when we consider large a, there are still allowed 
regions which can be probed by Advanced LIGO with both burst detection and the GW background. 
Choosing parameter values from such a region, we calculate the Fisher information matrix to study 
how the cosmic string parameters can be constrained by future GW experiments in the following 
sections. 



4.3 Formalism of Fisher analysis 

The maximum likelihood method is widely used to estimate model parameters in the analysis of 
cosmological observations [58]. For a given data-set, the set of parameters that is most likely to 
result in the model prediction are those which maximize likelihood function C The error in this 
estimation can be predicted by calculating the Fisher information matrix, which is defined as 

_ dHnC 

Assuming a Gaussian likelihood, the expected error in the parameter 6i is given by 

ere, = Vi^i (40) 

Here, we apply this Fisher matrix formalism to predict how accurately cosmic string parameters 
can be constrained in future direct detection experiments, for both cases of burst and stochastic 
background detection. 
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4.3.1 Burst detection with a single detector 

If GW bursts from cosmic strings are detected frequently by GW detectors, we would be able to 
make a catalogue of the bursts from cosmic strings. Let us suppose that we have a large enough 
number of samples and the number of observed GW bursts per strain interval hi to hi + dhi is given 
by 

Ni = ^{K)dhi, (41) 

where ^{h) = dR/dh x T is predictable from the cosmic string parameters as presented in Sec [31 
Depending on the detector sensitivity, the catalogue has a magnitude limit oi h = hmm- 

Assuming that the number of GW bursts follows a Poisson distribution [501 EHl EOl [61], the 
probability of observing k events in each strain bin is given by 

P. - (42) 
The likelihood function is defined by the total probability of all bins as 

i 

Substituting the likelihood into Eq. (!39|) . we obtain the Fisher matrix 

^de^de^N,' ^^^^ 

where we have used k Ni. Substituting Eq. (HTl) and rewriting in terms of integral, the Fisher 
matrix is given by 

p a<j) a$ 1 

We use bursts detectable by the future GW detector network with p > 4:, which corresponds to 
the limit of the catalogue being //imin 

~ 3.4 X 10"23 at / = 220Hz, as mentioned in Sec. liXTl |^ 

4.3.2 Search for the stochastic GW background 

In the case of cross-correlation analysis, assuming the Gaussian distribution of the data /t around 
the mean value (/i), the probability distribution function (likelihood) is given by 



C = , exp 



2^2 



(46) 



^ Higher signal to noise ratio may be required for the use of the Fisher matrix approximation, since the distribution 
of the ampHtude A deviates from the Gaussian shape because of uncertainties in determining the sky location of the 
burst [54] ■ Here, however, we assume that the sky locations of bursts are determined with a good accuracy by using 
the multiple detector network. 
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From Eq. f l5^ . we obtain the Fisher matrix 

a2 d9i ddm' ^ ^ 

Substituting Eqs. fl32p and f l36p . and the optimal function used in Sec. 14.1.2^ the Fisher matrix (for 
detectors between I and J) is rewritten as 

•^-'--^lO^J 2Ty^ df ;a5„^^(;)5„^^(;) ' (48) 

For multiple detectors, the Fisher matrix can be calculated with 

3HI y a; p \lij{f)\'de,nGw{f)de„^^Gw{f) 



■^"»=iio^j '^^^X rs„M)s..Af) • 

4.4 Constraints on cosmic string parameters 

Using the Fisher matrix formalism presented above, we forecast constraints on cosmic string param- 
eters from future direct detection experiments. In Figure El we present an example of the expected 
future constraints in the case where the parameters are Gfi = 10~^, a = 10~^^, p = 1. Each ellipse 
represents the 2a error contours expected from 3 years of observation with future ground-based GW 
network (Advanced LIGO, LCGT, VIRGO, and AIGO with the sensitivity given in Eq. ([25])). In 
this setup, 62724 bursts are detected with p > 4 by multiple detectors in 3-year observation, and the 
GW background is detected with p ~ 57 by correlating outputs of all the detectors. We clearly see 
the constraints from burst detection (black line) are tighten when it is combined with that from the 
GW background (red line). 

The reason why we can obtain the stronger constraints by combining burst detection and GW 
background measurements is that the constraints from the two measurements have different param- 
eter degeneracies. In the case of the burst detection, the direction of the parameter degeneracy is 
determined by the parameter dependence of the rate dR/dlnh, which is presented in Appendix Rl 
The parameter set used here corresponds to the case (iii), i.e. Eq. (153|) . with //i3,3 ~ 4.7 x 10~^^. 
Since GWs larger than ^13^3 are observed via burst detection, the direction of the parameter de- 
generacy is oc {Gjjya^^^p~^. In the case of the GW background, the direction of the parameter 
degeneracy is determined by the parameter dependence of flew, which is investigated in Appendix 
iBl In this parameter set (the case (iii)), the GW background is given by (160 p . and the second term 
dominates the first term if it is evaluated at the frequency of Advanced LIGO. Thus, the constraints 
from the GW background has the parameter degeneracy of oc Gfia~^^^p~^. The constraints from the 
GW background have strong parameter degeneracies since the observable is basically only Qgw at 
/ = 220Hz. By combining it with the constraints from bursts, the degeneracies are broken. 

Here, we would like to mention when these bursts and the background detectable by Advanced 
LIGO are emitted in the case of the parameter set we take above. As mentioned in section 13.21 



^ Note that, the direction of the degeneracy seen in the figures does not directly corresponds to the parameter 
dependence described here, since the shown constraints are marginahzed over the other parameter. 
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Figure 6: Marginalized 2a constraints on cosmic string parameters in p — a, a — G/i, and p — Gfi 
planes, respectively. The fiducial parameters are taken to be G/x = 10"'', a = 10~^^, p = 1. The 
solid black line represents the constraints from the burst detection alone, and the red line represents 
the combined constraints from the burst detection and the GW background measurement. 
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and Appendix [XJ bursts of given frequency / and amplitude h are mainly emitted at redshift which 
satisfies l{f,h,z) ~ at{z) for a < TGfi, so bursts detectable by Advanced LIGO are generated 
at 2; <^ 1. To be more specific, GWs of amplitude equal to the sensitivity of Advanced LIGO, 
fh ~ 3.4 X 10"^^, and of frequency / = 220Hz come from z ~ 3.6 x 10"^, and those detected 
once per year by Advanced LIGO, which have fh ~ 10~^^ and / = 220Hz, come from z ~ 10""^. 
The main contribution on the GW background comes from GWs emitted at redshift which satisfies 
/(I + 2) ~ at{z) in the radiation-dominated epoch and those emitted at z ~ 1. At the frequency of 
Advanced LIGO, the contribution from the latter is larger. 

5 Conclusion 

Future GW experiments can be a unique and useful tool to test the existence of cosmic strings. If 
GWs from cosmic strings are detected, they could provide important constraints on cosmic string 
parameters. In this paper, we have studied the potential of upcoming ground-based GW experiments 
to search for signals from cosmic strings and estimated the power to place constraints on cosmic string 
parameters. 

The key point of this paper is that GW experiments can search for cosmic string signals both 
in the way of burst detection and GW background measurements. First, we find the parameter 
region where GWs from cosmic strings are detectable by future experiments considering the both 
cases. Furthermore, we investigate constraints on cosmic string parameters obtainable from direct 
detection measurements of the GW bursts and the GW background, and found that their information 
is complementary from each other. Thus, if both the GW bursts and the GW background are 
detected, we can tighten the constraints on cosmic string parameters by combining data from the 
two different measurements. Although we demonstrate only the case of ground-based experiments, 
this is also the case in future satellite experiments like BBO and DECIGO which is designed to search 
for both GW bursts and GW backgrounds. 

One thing we must note is that there would be another sources of GW bursts and GW back- 
grounds. There are surely many astrophysical candidates which generate GW bursts. However, since 
GW bursts from cosmic strings have a characteristic frequency dependence, cosmic string bursts 
could be distinguished from those from other sources. In the case of the GW background, although 
there is no certain source, some models can predict a GW background around the frequency band of 
LIGO sensitivity. We could also use information on the frequency dependence of the spectrum, but 
it may be difficult to identify whether the detected GW background originates from cosmic strings 
or other models. However, although it depends on the values of the cosmic string parameters, we 
may be able to use other observations which explore different frequency ranges of GWs such as CMB 
B-mode measurements [65l EH [671 EH EH [TQl [7T1 [72] and pulsar timing experiments [551 U3\ to verify 
the GW background from cosmic string, which will be studied in our future work. 
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Appendix 



A Parametric dependence of the burst rate 

In this appendix, we roughly estimate the dependence of the burst rate on cosmic string parameters. 
Cosmic strings expiring at each epoch give dominant contribution to GWs, because they are the 
most numerous. So, roughly speaking, we only need to estimate their contribution to the burst rate, 
as shown formally below. Since we are interested in high frequency GWs detectable by ground-based 
GW experiments, we consider bursts which satisfy / > (minja, FG/i})"^ t^^(2eg)(l + z^q)"^, which 
means that expiring loops start to contribute to relevant frequencies in the radiation-dominated era 
and their contribution continues until today. 

The rate dR/dlnh is derived by integrating Eq. f|T6l) in terms of z. Using Eqs. f|TT]) . f fT7|) . and 
(|T8|) . we find that in the case where a ^ TGfi, so loops are long-lived, the dominant contribution to 
the integration of Eq. fITB]) comes from the redshift around Zm which satisfies 



for large h. The size of loops which emit GW bursts of frequency / and amplitude h depends 
on when they are emitted and it is denoted by l{f,h,z). Equation ( l50i) means that, for a given 
frequency and amplitude, GWs are mostly generated by loops expiring at redshift Zm- Therefore, 
we can roughly estimate the rate by dR/dh ~ Zmif,h) x -^^^{f,h,Zm{f,h)), neglecting accuracy 
of 0{1) numerical factors. Note that for too small h, there is no solution which satisfies both 
Eq. (150!) and /(I + z^) > l^^{f,h,Zm), which corresponds to the lower frequency cutoff of bursts. 
In this case, the main contribution to the z integral of Eq. (1161) comes from Zm which satisfies 



In the case of short-lived loops, a <^ TGfi, Zm is given by /(/, h, Zm) = at{zm). In this case, if we 
consider small h, any value of z does not satisfy /(/, h, z) = at{z), so dR/dln h is strongly suppressed. 

The rate of GW burst depends on when bursts are mainly emitted and when loops which emit 
such bursts are formed. We summarize the results below: 



(Here, Qr is the current energy density of radiation divided by the critical density.) 

In this case, loops are so long-lived that those formed in the matter-dominated era does not 
decay before today. In other words, loops which have decayed by today are formed in the radiation- 
dominated era. 

The burst rate is given as follows. 

^ We also neglect the dependence on F in the following estimation, since it is expected to be canceled by numerical 
factors which determine F, for example, that in Eq. (|12p . 



/(/, h,Zm) = TGfit{Zm), 



(50) 



f{l + zj = l-\f,h,zj. 
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dR 
dlnh 



if,h) 



(G/.)l/16al/2^-2 (^Y^'^ /-13/4i-37/16^^_31/16 

where to is the present age of the Universe. 

Bursts whose amphtudes are hi i < h < hi 2, /ii,2 < h < /ii 3, and h > /ii 3 correspond to those 
emitted in the radiation-dominated era, in the matter-dominated era, and at redshift smaller than 1, 
respectively, li h < hi^i, any Zm does not satisfy /(I + Zm) ~ ^~^{f, h, Zm)- Therefore, bursts which 
have amplitude of h < hi i are emitted in the radiation-dominated era not by loops expiring but by 
loops which still have their lifetime. 

Increase of Gn enhances the amplitude of each burst. Since bursts of a given amplitude can come 
from earlier epochs and more distant points for larger Gfi, their number increases. This is why the 
rate is proportional to a positive power of Gfi. However, there is also an effect which decreases the 
rate with increasing Gfi, that is, reduction of the loop lifetime and the density of expiring loops due 
to the enhancement of the efficiency of the GW emission. Note that loops formed earlier are more 
abundant than those formed later despite more dilution by the cosmic expansion, since the density at 
their formation is larger for older loops. The total power of Gfi is determined by these effects. As a 
increases, the elongation of loop lifetime enhances the rate but the initial number of loops decreases. 
In this case, the former is more efficient and the power of a in Eq. flSTl) becomes positive. 



-1/2 



Efl = "-1,1 



11/5 



(51) 



(ii) TGfi <a< TGfi (^^) 



3/2 



In this case, loops which expire at small redshift are formed in the matter-dominated era. The 
burst rate is given by 



dR 
dlnh 



if,h) 



-1/4 



(G/.)V2aV2^-2(r^) ' fV2t^k-i/2 ;/,<(^) ' f-H^ 



(G/i)V2al/2^^-2 (^^^ ^^/^°^-18/5^-12/5^_ii/5 

■,h,,^ < h < {G^,r/'f (^)"'^' = h,,, 

(G/i)Vl6al/2^-2 (^^)"'^' j-13/4^-37/16^_31/16 

;/i2,i <h< {Gf,y/^a^/^r^/h-'/' {wY"^' 
{Gf,y/^^~'r^'/%'/'h~^ ■,h>h2,. 



-1/2 



h 



2,1 



(52) 



'-2,3 
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Bursts whose amplitudes are /i2,i < h < /i2,2 and /i2,2 < h < /i2,3 correspond to those emitted 
in the radiation-dominated era and in the matter-dominated era by loops formed in the radiation- 
dominated era. On the other hand, bursts in the range of /i2,3 < h < /i2,4 and h > /i2,4 are emitted 
at redshift 1 < z < Zeq and 2; < 1 by loops formed after the matter-radiation equality. These are 
mainly emitted by expiring loops. Bursts of h < /i2,i, on the other hand, are emitted by loops which 
do not expire soon. 



(iii) a < VGji 

In this case, loops are short-lived. The rate is given by 



■,G^.a-'rhf [^y"" ^ /.3,i < h < G/zaV3/-4/3,-V3 ^n^y'/' ^ ^^^^ 

(G/.)2al/3^^2^-14/3i-8/3^„3 .^^j^^^^ 

Bursts in the range of h^^i < h < h^^2i ^3,2 < h < h^^^, and h > h^^^ are emitted in the radiation- 
dominated era, in the matter-dominated era, and at z ^ 1, respectively. The burst rate for h < h^ i 
is suppressed for the same reason as /i < hi i and h < /i2,i. In these cases, elongation or shortening 
of loop lifetime do not affect the rate and the rate depends only on parameters through time of burst 
emission and initial number of loops. 



dR 
dlnh 



B Parametric dependence of the spectrum of the stochastic 
GW background 

In this appendix, we roughly describe the amplitude of ^qw in terms of cosmic string parameters. 
In the same way as the previous section, we consider only high frequency GWs. 

(i) a > TGfi {%^) 

In this cases, we find through Eq. (15T!) that the dominant contribution to the integral of h"^^ in 
Eq. (12T|) comes from h ~ hi^i or /i ~ /113. Approximating that the contribution to J dlnhh^^ only 
comes from these GWs simply as hf^^{hm), where hm is hi i or /ii 3, we get 

Ogw - (Gf^y^'^'^V^ + (G/.) V^aV (^) t^'^'f-'^'. (54) 

The first term, the contribution from h ~ hi^i, represents the GWs emitted when loops decay in the 
radiation-dominated era, and the second one, the contribution from h ~ /ii,3, corresponds to GWs 
emitted at 2; ~ 1 by loops which expired recently or are decaying today. 
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We can estimate that the first term is equal to the energy of loops expiring at time t (redshift z) 
which satisfies TGfit ~ + z) as 



^ Ploopit)\l^rGfit 

i^GW 7TT 7^ 

Ptot (, t j i 

Ploopi^i^ll^ati ^ ~\~ Zi 

r\_/ 

PtotiU) 1 + Z^m 

/ g \ Qr_ 



\TGfI J ilr^ 



~ iG^^y/'a'/^^;'^. (55) 

Here, ti and Zi are the time and redshift of the loop formation respectively, which is related to 
the time of GW emission as t ~ ^^t, and the factor ^ denote the dilution of GWs from the 

* a ' Urn 

matter-radiation equality to now. 

The second term can be estimated as follows, 

Ogw - ^(1 + z,r X ^(/)/ X to X ^ ~ (G/i)VV/V' ^§^V'^^o'^V"^/^ (56) 

Trtt^i df Per J 

where ti ~ ^^^o and Zi are the time and the redshift when loops expiring today are formed. Note 
that the number of GWs emitted recently is so small that it might not be regarded as a part of the 
background. The condition -^y^{h = hi^s) > f leads to 

(G/.)"i^/^aV (^) if to)-'/' > 1, (57) 
which can be satisfied by small G/i, for example Gfi < 10^^\ for / ~ lOOHz, a = 0.1, p = 1. 

/ \ 3/2 

(ii) TGfi <a< TGfi f 

One can get ficw by estimating J dhh"^^ in the same way as demonstrated above. The result is 

^^Gw - (G/.)^/^a^/V^#^ + iG^,f'i;n%'^'r'/'- (58) 

Again, the first term corresponds to the GWs emitted by loops decayed in the radiation-dominated 
era and the second one is the contribution from loops expiring today. The first term is identical to 
the one in Eq. ( 154|) . On the other hand, the second one is different because loops expiring today are 
formed after matter-radiation equality in this case. 
The condition to include the second term in flaw is 

iGp)-'/'l-Vtor'/' > 1- (59) 
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For / ~ lOOHz and p = 1, this reads Gfi < 10 



(iii) a < VGji 



In this case, we get 



1/3^-1/3 r-l/3 



(60) 







Again, the first term is the contribution from GWs emitted in the radiation-dominated era and the 
second one represents GWs emitted recently. The first term coincides with the ratio of the energy 
of the infinite string network to the total energy, ~ GyU7~^, except the dilution factor because, 
in this case, loops released from infinite strings immediately expire by radiating GWs. The second 
term also consists of GfX'y~^ and the factor {fato)~^^^, which corresponds to the tilt of the spectrum 
of emitted GWs. 

The condition to include the second term in Qqw is 



which is relatively easy to satisfy for small a. For example, this reads Gfi ^ 10 ^ for / ~ lOOHz, 
a = 10~^^ and p = 1. 
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